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We discuss modifications of the thermal Dark Matter (DM) relic abundances in stringy cosmolo- 
gies with D-particle space-time foamy backgrounds. As a result of back-reaction of massive DM on 
the background space-time, owing to its interaction with D-particle defects in the foam, quantum 
fluctuations are induced in the space-time metric. We demonstrate that these lead to the presence 
of extra source terms in the Boltzmann equation used to determine the thermal dark matter relic 
abundances. The source terms are determined by the specific form of the induced metric deforma- 
tions; the latter depend on the momentum transfer of the DM particle during its interactions with 
C^l . the D-particle defects and so are akin to Finsler metrics. In the case of low string scales arising from 

large extra dimensions our results may have phenomenological implications for the search of viable 
jjj '" supersymmetric models. 
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I. INTRODUCTION 



The nature of the Dark sector of our Universe constitutes one of the major unresolved puzzles of modern physics. 
Indeed, according to observations over the past twelve years, 96 % of our Universe energy budget consists of unknown 
entities: 23% is Dark Matter (DM) and 73 % is dark energy (DE), a mysterious form of ground state energy. DE 
is believed to be responsible for the current-era acceleration of the Universe. These numbers have been obtained by 
best- fit analyses of a plethora of astrophysical data to the so-called Standard Cosmological Model (ACDM), which 
is a Friedman-Robertson- Walker (FRW) cosmology, involving cold DM, as the dominant DM species, and a positive 
cosmological constant A > 0; the data range from direct observations of the Universe acceleration, using type-la 
supernovae [lj, to cosmic microwave background @, baryon oscillation [3| and weak lensing data Q. It should be 
stressed that the afore-mentioned energy budget depends crucially on the theoretical model for the Universe considered. 

An interesting, and not commonly discussed, class of Cosmological models that may lead to modifications of the 
Dark sector, involves space-time with a "foamy" structure at microscopic (Planckian or string) scales [5j, due to 
quantum gravitational interactions. In the past, for a variety of reasons, such models (differing in the details of the 
constructs of space-time foam) have been considered by many authors. They exhibit a profusion of features that 
can be falsified experimentally and their predictions range from light-cone fluctuations caused by stochastic metric 
fluctuations @, to macroscopic Lorentz symmetry violations 0|- In this note we would like to present a first study 
towards the contribution of such space-time metric stochastic fluctuations on the Dark sector and in particular on 
the DM sector of the Universe. We shall focus on a particular class of stochastic space-time foam models, inspired 
by certain types of string theory. These involve localized space-time defects (DO-branes or D-particles) [8] which are 
either allowed background configurations [8, 9] or arise effectively from suitable compactifications of higher dimensional 



branes (e.g. D3-branes wrapped up in appropriate three cycles in the context of type IIB strings (10|) . Observers 
on the brane detect a foam-like structure due to the crossing of the brane by D-particles. In this higher-dimensional 
geometry, only gravitational fields are allowed to propagate in the bulk; all other particle excitations, including DM 
candidates, are assumed to be described by open strings with their ends attached to the brane world. The brane is 
assumed to have three large spatial dimensions, and - depending on the model of string theory considered - it may 
have a number of compactified extra dimensions. 

Dynamical D-particles should not be viewed as material excitations of the vacuum but rather as vacuum struc- 
tures. This contrasts with attempts to represent such D-particles as ordinary superheavy DM excitations from the 
vacuum (Tf| , owing to the completely localized nature of the D-particles in the extra dimensions. In our construc- 
tion, closest in spirit to weak coupling string theory, they are just vacuum defects. In fact it can be shown that the 
gravitational interactions among such (BPS) D-particles are cancelled by appropriate gauged repulsive forces induced 
on them from other branes in our supersymmetric models of D-foam [8|. Hence such a collection of D-particles in 
the bulk does not affect the Hubble expansion on the D3-brane worlds, and so their concentration cannot be re- 
stricted by considerations on overclosure of the Universe within our models. Important constraints on the density of 
defects in D-foam models can still be imposed by astrophysical experiments on the arrival times of high energy cosmic 
photons According to these string- foam models, there should be an effective refractive index in vacuo, such 

that higher energy photons would be delayed more, since they would cause stronger disturbance on the background 
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space-time The interaction of material open strings (such as photons) with the D-particles leads to an induced 

distortion of space-time described by a metric, which depends on both the coordinate and momentum transfer of the 
photon during its scattering with the defect and so has similarities to a Finsler metric [l3j]. This is a topologically 
non-trivial process, involving the creation of a non-local intermediate string state, oscillating from zero length to a 
maximum one, according to a time-space stringy uncertainty [^.[l4|. This causes a time delay for the photon emerging 
after capture by the defect, proportional to the incident energy of the photon. . 

The purpose here is to analyze, in the same spirit, the modification of the estimate of the DM budget of the Universe 
as compared to the ACDM model due to the quantum fluctuations of the D-particles. In fact, as we shall argue, the 
propagation of massive DM particles on such space-times, induces a back-reaction, which in turn has consequences 
for the amount of thermal DM relics of these particles; this, in turn, impacts on astro-particle tests of particle 
physics models incorporating supersymmetry, that provide currently one of the leading DM candidate species, the 
neutralino. Its thermal abundance, calculated within the simplest supersymmetry models (minimal supersymmetric 
model embedded in minimal supergravity (Isj). is heavily restricted by cosmic microwave background data; hence the 
available parameter space for these simplest supersymmetric models may vanish in the near future, on incorporating 
also data from collider experiments such as the LHC at CERN [l^. These constraints depend strongly on details 
of theoretical models. In the presence, for instance, of extended supergravity models with time-dependent dilaton-</> 
sources [ItJ, the calculated amount of thermal neutralino relic abundance can be smaller than the one calculated 
within the ACDM-minimal supergravity cosmology. Such dilaton models allow more scope for supersymmetry, which 
can thus survive the otherwise stringent tests at the LHC |18j. In our work we will find that the effects of the D-foam 
on the thermal relic abundances oppose those from the dilaton models. These effects may become relevant for models 
allowing low string mass scales. 

For clarity we will first outline the basic reason behind such modifications in the DM thermal abundances. Non- 
equilibrium cosmology models are associated with space-time distortions, due to either the presence of time dependent 
dilaton sources (cf supercritical (SSC) dilaton quintessence string cosmologies pj]]), or the induced back reaction of 
the DM particles onto the space-time itself. Boltzmann equations are used to determines the thermal DM species 
cosmic abundances. The effect of space-time distortions can be subsumed as extra contributions to a source T(t) 
for particle production (at cosmic Robertson- Walker time t) on the right-hand side of the appropriate Boltzmann 
equation. In general (for a Universe with three spatial dimensions (i.e. a 3-brane), with no loss of particles to the 
bulk ), with n(t) the density of the DM species, the apposite Boltzmann equation reads : 

d ti 

— + 3Hn = T{t)n + C[n] (1.1) 

da 

where H = — = - is the Hubble ratio, a being the scale factor in the Robertson- Walker metric, and C\n] is the 

a a ' ° i L J 

standard collision term describing the deviations of the species population from thermal equilibrium. In standard 
cosmologies, which we assume here, this has the form [19( 

C[n} = -(av)(n 2 ~n 2 cq ) (1.2) 

with n eq a thermal equilibrium density of species, a the total cross section evaluated in the background metric, and 
v the Mceller velocity. 

We will provide an explicit expression for the D-foam induced source T(t). The interaction of DM particles, 
represented as open string excitations on the D3 brane world, with D-particles, (interpreted as space-time defects in 
the ground-state ), results in metric distortions. The latter depend on the momentum transfer of the matter particles 
and so are of Finsler type [13j. It was observed that [20j, for the universe to expand by a scale factor a(t), in a way 
consistent with the world-sheet logarithmic conformal invariance of the D-particle recoil process during capture and 
scattering (2l| . there is an induced space-time distortion which has the form: 

ds 2 = -dt 2 + ma 2 (t )dtdx l + a 2 (t)dx i dx l S. lJ , t > t . (1.3) 

Here the underlying space-time is given by a spatially-flat Robertson- Walker metric, which is the space observed by 
a low-energy observer on the brane world. In (|1.3[) . g s is the string coupling, M s is the string state, and M s /g s is 
the mass scale of the quantum gravitational foam fluctuations in the model. The metric (|1.3[) has precisely the form 
corresponding to a boosted frame (the D-particle's rest frame), with the boost occurring suddenly at time t = to, 
the time of the capture of the string state by the D-particle. On account of momentum conservation during the 
scattering (2l| . the D-particle recoil velocity Vi is related to the momentum transfer Api of the string state, which in 
turn can be parametrized as a fraction of the incident moment um (22l| 

^Pi £iPi 100 ^1 /t a\ 

Vi = 9s-rr- = nPi = 9s-rj- , no sum over i = 1, 2, 3 , r, < 1 , (1.4) 

1 vl s 1VI s 
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where pi is a co-variant co-moving momentum of the DM particle in the (distorted) Robertson- Walker background. 
In our foam model we assume that the population of D-particles, is approximately uniform and relatively dense over a 
given epoch of the Universe. In such a case one may average (|1.3j) using appropriate distribution functions. Denoting 
such an average by <C • • • S>, and assuming stochastic Gaussian distributions, 

< n »= , < TiTj >= of % , no sum over i = 1, 2, 3 . (1.5) 

The variances af need not be independent of i; however to keep in line with the observed isotropy of the Universe 
at large scales, such potential anisotropies (due to an anisotropy of the population of D-particles in the bulk [gl, Q) 
would have to be small. 

Upon considering such foam populations, where the time scale between D-particle captures and emissions is much 
shorter than the inverse mass of the DM, a coarse-grained description is applicable; the propagation of a DM distri- 
bution over such a background is described on average by metrics <C ds 2 ^> in which to a good approximation to ~ t. 
Consequently we write the distorted metric in the following form: 



f),,,. 



/ -1 a 2 (i)ripi a{t) 2 r 2 p 2 a 2 (t)r 3 p 3 

a 2 (t)r lPl a 2 (t) 

a 2 (t)r 2 p 2 a 2 (t) 

\a 2 (t)r 3 p 3 a 2 {t) 



(1.6) 



The averaging over foam populations is performed using (| 1 . 5j) ; hence, when considering the evolution of DM densities 
in such backgrounds, terms with an odd number of r^'s will be ignored. This should be understood in our subsequent 
discussion. The momentum dependence of the metric (|1.6|) implies a modification of the form of the geodesic equation 
for ^ppr together with the modification of the Christoffel symbols. Keeping terms up to order r 2 and dropping 
everywhere cross terms of the form r^j for i =/= j (these terms will in any case vanish when one takes the average), it 
is straightforward to show that the pertinent geodesic equation for /j, = i (= 1, 2, 3) reads: 

+ ^(ww (p ) 3 - ^mrbVY. (^') 2 (1-7) 

j j 

where the Einstein summation convention is not applied. 

To discuss DM relics we need to solve the relevant Boltzmann equation that describes the evolution of the phase- 
space density function of the DM species in t. The distribution function / of a particle species is specified by phase 
space variables of the system; it is thus natural to define a local space-time momentum in an expanding universe (23| : 

p l = a{t)p\ i = 1,2,3 . (1.8) 

In terms of these scaled momenta, the energy-momentum dispersion relation for a DM particle of mass to in 
our (spatially flat Robertson- Walker space-time background), assumes an effectively "Minkowski-space-time" form, 
p^p v g^ v = —E 2 + a 2 (t)p z p : '5ij = —E 2 -\-p l %P8ij — —to 2 . Thus, it is essential to define the phase-space densities as 
functions of the coordinates z M = (x\t) and the local momenta p\ f(x fi ,p l ). This choice of variables is particularly 
important in keeping the correct scaling properties of the DM density with the scale factor. In the usual isotropic 
Robertson Walker background, the momenta are assumed on-shell, and so the phase-space density depends on|p| (i.e. 
E) rather than p. However, in our case, the small anisotropies that characterize the foam fluctuations are taken into 
account by assuming a dependence on the individual components p { . 
The Liouville operator acting on / (x^, p 1 ) takes the form: 

L[/]=/^ +m y^f. (i.9) 

i 

We apply the isotropy condition , ^ = 0, i = 1,2,3, and note that ^ = a(t)^r + d(t)^p l = a(t)^- + a{t)^p l ; 
the following expression for the Liouville evolution operator is obtained: 

^ = ^--HYj?—- 2Ha 2 (t)p° Y n f *L + 8Ha 4 (t) V r 2 (pf °l 
j3° dt ^ y dp 1 K 1 ^ df K ' *f V ' dp 1 

j i i j i 
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The mass shell condition gives the energy dispersion relation 



P = a 2 (t)j2(pfr i+ lj2(fy 



« 4 W (Eiiffn)' 

1 + — - 



(1.11) 



Further analysis involves the approximation of heavy DM, which we consider in this work, as the phenomenologically 
dominant species are the heavy ones that in general are expected to cause the most significant distortions in the space 
time background: m 2 3> X^P 1 2 ■ ^ n this case expanding (|1.11[) up to second order and averaging over the ensembles 1 , 
for the random variables fj, we obtain for the Boltzmann equation (now with the binary collision term taken into 
account): 

I - H E + E *? § - E E (^) 3 § + 4ffmv ( t) E 

-4^5:^) 2 E^¥f (i.i2) 

3 i 

To make the connection with 11.11 the number density n (t) is defined as 

n{t)^^ J d 3 pf(trf) (1.13) 

where e? 3 p = dp 1 dp 2 dp 3 . Following [23| an average temperature T is defined through the equation 

9 ' d 3 p (p 1 ) 2 f= Tmn. (1.14) 



(2tt) 3 



In the case of superheavy dark matter, the fourth and sixth terms in (|1.12j) are small compared to the others (and 
so will be neglected). Integrating (|1.12l) with respect to d 3 p, yields: 

£ + 3*„ = r(«,„ + ^/^f! (1.15, 



where we have incorporated all extra corrective terms in a time-dependent source term given by: 

T(t) = Ha 4 (t) (E ^ [ 1 8Tto + 4to 2 ] , (1.16) 

and 



,3 £[/] 



(2tt) 3 J r E 

The reader should bear in mind that in our analysis we have kept only leading order terms in the fluctuations. 
Moreover, the correct scaling of n with the Hubble parameter, embodied in the term 3Hn on the left-hand-side 
of (|1.15p . has been obtained because the local momentum p> has been used in the argument of the phase-space 
distribution function f(x l ,t, p J ) [23], as mentioned earlier. 

Hence the role of the stochastically fluctuating induced Finsler metrics in the D-particle foam model is simply to 
give rise to particle-production source terms in the Boltzmann equations, linear in the DM density. The presence of 
these terms (or equivalently the metric distortions) will affect quantities such as the thermal relic abundances (which 
will be distinguished with a prime). Although formally the situation is similar to that of time-dependent dilaton and 
non-critical string sources encountered in the super critical string (SSC) cosmology model of [17] . our modification 



The reason why we can do this at the equation level is because the time scale of D-particle scatterings is assumed to be much shorter 
than the Hubble time H . 
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£l' x of the thermal relic abundance of a single heavy DM species \ wm carry extra terms as compared to [TtJ ■ The 
process we follow is presented analytically in [25]. In a standard notation, the result is: 

n' x h 2 f x ° T{x) ^ 1 




= 1+ -Wzdx-- J(y) ^dx dy (1-17) 



(^X^o)no source \ J Xf Hx J J Xf \JlVf ^ X J J I \9eff J x=x . 



where J = J (va) dx, J(x) = ^ v %' and x = m x /T, following standard usage; x$ denotes the current value of x 

Xq 

corresponding to the CMB temperature Tomb — 2.7° K, and Xf denotes the (D-foam-modified) freeze-out temperature, 
estimated by using the freeze-out criterion Y (xf) — Y cq (xf) ~ coY cq (xf), where Y(x) = 2152, with s an entropy 
density (satisfying 4? (sa 3 ) = 0) and the suffix "eg" denotes the equilibrium expressions (in the presence of the 
deformed metric) whose analytic form is given in (25j . The quantity Co is a phenomenological constant of order O (1), 
which can be determined by appropriate numerical fits. An analytic expression for Xf is difficult to obtain in our 
case [25j; however, for weak foam effects, approximate expressions can be found by an appropriate expansion 2 . The 
quantity (fi x /io)no SOU rce = 1,066x] A f GcV — ^jQeffJ denotes the relic density in the absence of a source term T = 

Xf 

and J = J (va) dx. However, the reader should bear in mind that what we denoted as (f2 x /ig) no SOU rcc is not quite 

xo 

the standard expression within the FRW Cosmology, since (va) carries implicit information about the D-particles 
effects on the total cross section (denoted by a prime) and also the freeze-out point that appears as the upper end of 
integration is shifted as described above; Mp is the Planck mass, which is related to the D-particle mass scale M s /g s 
via 

M 2+s 

M P = ^WV} S \ (1.18) 
9s 

with Vc^ the volume of the compactified <5-extra dimensions of the 3-brane world, which depends on the specific 
model under consideration; the Mp (= 1/Gjv) enters the formalism through the effective four-dimensional Einstein- 
Friedmann equation H 2 = 87r ^ w (p + Ap). Here, p denotes the total (critical) radiation and matter energy density, 
including DM, and Ap symbolizes collectively all the stochastic effects of the D-foam on the standard Einstein's 
equations. Regarding the latter, we note a small but important difference compared with the case of the time- 
dependent dilaton model of [17j. The stochastic fluctuations of the space-time due to the recoil- velocity fluctuations 
of the D-particle defects in the D-foam, do contribute to thermalization of the Universe, indirectly, as a result of 
their coupling with the photon or electrically-neutral matter excitations (such as DM) via the distorted Finsler metric 
backgrounds 3 . 

In order to calculate the explicit form of the factor ( j that contributes to the D-foam corrections to dark 

\9eff J x — x ^ 

matter relic abundances, one should first derive the equilibrium expressions for relativistic fermions and bosons in the 



2 For instance, it can be shown that the freeze-out criterion can be written as f25j : ex P (/J^ ~ ' dx'j -g^p-dy = (cq + 
l)~ 1 uoa;J 3 n^j 1 (x f) where T(x) is the source 1 11,161 . the total cross section cr(x) ~ 5qX~ j = 0(1) for s(p)-wave annihilators [lSl . 

uq = ^g-h' m 3 , where h' denotes the entropy degrees of freedom, and we used H = H m x~ 2 , H m = X.&Tg'^J^ , with Mp the Planck 
mass; x in corresponds to some initial value of the (inverse) temperature, say at the end of the inflationary era. The equilibrium number 
density n cq (x) receives corrections from the foam |25|| . For weak- foam situations, corresponding to large values of M s /g s compared 
to the DM mass m, one may expand the exponential to first order in T(x), taking into account that the scale factor a(t) ~ 1/T (c.f. 
below) , and solve iteratively for x j . On assuming that the fitting constant cq can be chosen in such a way that the foam corrections on 
the right-hand-side of the above freeze-out criterion equation are subleading compared to their counterparts on the left-hand-side, one 
finds self-consistently |25|| that the freeze-out point Xf slightly increases relative to its foam-free-case value, x^ (which for neutralino 

dark matter is of order 20 [H): x f j) ~ x ( f 0) ' U) + (g„ J I J I ^— ) fef =1 A 2 ) , where j = (1) for s- (p-) annihi- 
lators, and A? is a dimensionless foam-fluctuation variance [c.f. 11.41 1). Hence, the foam leads to a relative decrease of the freeze-out 
temperature. However, taking into account that for neutralino masses in the range of O(10 2 ) GeV [H, x in ~ 10~ 12 (since one can 
reasonably place the end of inflation at temperatures of order 10 14 GeV), and today's CMB value xq ~ 10 14 , one observes that the shift 
in the freeze-out temperature is neg ligible, for A 2 < 1, even for low string scales of order TeV. 

3 In our D-foam models [8l 4lOH charged matter cannot couple dominantly to the D-particles due to electric charge flux conservation. 
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presence of the D-foam. To this end, we use the basic formula |25| : 

P = 7A3 j «™ r » d 3 p (1.19) 

where uj t is the equivalent of (jl.lip with a(t) = 1 and the notation <C . . . ^> denotes the average with respect to the 
statistical parameters of the foam i.e. : 



1 r r ? 

< ncj r >= TT = / drj < nuj > r cxp(--4r) (1.20) 

. cr 7 - V27T 7 ^ 



^ —00 

and (for any member of the ensemble of possible foams) we have the canonical energy distribution 



< nu> > r = 



exp(/?(w r — (j,)) + zd 



where w = +1 (—1) applies to fermions (bosons) and j3 = -j-^ . 
In this way we then obtain the expression (25l j: 

" V 4 ^/ 3 (l,/3M,^)+^ 2 r 6 ^/5(l,j8A* J w)-^/6(l,^,ro) + ^/ 6 (2, j 8/i,w)-^/ 6 (l, J 8/i,tu) 

+ ^£/ 6 (2, fin, w) - ~/ 7 (l, /3/x, + ~£/ 7 (2, ro) + |/ 7 (1, ^, w) - |e/7(2, /9/*. + ^, w) H2I) 

where <? denotes the spin degrees of freedom of either the fermions or the bosons under consideration and we have set 
: cf 2 = <j\ + a\ + a\ . The functions fj appearing in (|1.21j) are defined through integrals of the general form: 

f°° k j 

f s {l,H,Z) = dk— — (1.22) 

Jo (e fc -^ + 

that can be related to £ and Dirichlet r\ functions. We also note here that since we are interested in relativistic species 
(contributing to s) we can set the chemical potential to (fi <C w r ) m t) 1 ■ 2 lj) . g' e jf, the modified effective number 

of degrees of freedom, is determined from p tot , the total relativistic energy density, since p tot = ^g'^fT^; one finds 
that: 

793.92 ^ / Tj s 




where g e ff stands for the standard effective number of degrees of freedom in the absence of D-foam corrections: 

i j 

The numerical factors appearing in (|1.23[) are the result of calculating the values of the integrals (|1.22[) . A detailed 
analysis of all these will appear in [25] . At temperatures of interest (since a freeze-out temperature is of the order of a 
few GeV for typical dark matter candidates with masses in the range m w lGeV — 10 4 GeV (26[) all the species in the 
standard model behave as relativistic matter. A counting of degrees of freedom, then, yields: gf — g qU arks + gieptons — 
6x(2 x 2 x 3) + 3x(2 x 2)+3x2 = 90 and g b = g g iuons+gEW+g v hoton+gmggs = 8x2 + 3x3+2+1 = 28. Substituting 
these values in equations (jl.23l) and (|1.24[) and using the approximations m k 1 we find: g e // = 106.75 and 
g' e ff = 106.75 + 22138(7 2 T 2 . These considerations imply that the D-foam corrections to the effective degrees of freedom 

will be of order |fi w 1 + 207.38ct 2 T 2 . The dimensioned variances af are naturally suppressed by the (square) 

of the heavy D-particle mass scale (cf. (|1.4p . with <C £ 2 3>= A 2 , i = 1, 2, 3 a dimensionless variance, that could be 
naturally up to 0(1)). Thus, we obtain: 

~™<M§4 (L25) 
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To determine the final corrections to the relic abundances we may expand the relevant expressions up to order erf. 
Expressing the result in terms of A?, we finally obtain 4 : 



n' x h% 



h 2 ) 

X'^O/ 110 source 



1 + 207.38^-^ 



1/2 



(1.26) 



2 ml 



fei=i A f ) > °- As noted in it is possible to 



where clearly the dominant correction terms are of order 
choose 

From the positive signature of the foam-induced corrections we observe that the relic abundances are larger than the 
corresponding ones evaluated within a standard cosmological model. In this sense, the space-time foam background 
will leave less freedom for supersymmetry at colliders [l7|. This effect is opposite to those induced by a time-dependent 
dilaton in supercritical string cosmologies [l7|. In the latter case, the negative signature of the pertinent corrections, 
implied that there was less DM available today as compared to the standard cosmology calculations; this could then 
lead to much heavier supersymmetric partners produced at colliders, such as the LHC, with falsifiable signatures [l8j . 
However, the corrections due to D-foam are in general small, as expected, and they can only be significant for low 
string scale models and heavy DM candidates. Indeed, the string scale M s is a free parameter in the modern version 
of string theory (cf Ql.lgp ). For traditionally high string scales (M s > O(10 16 ) GeV), in order for the D-foam effects to 
be significant one needs superheavy DM, with masses higher than M s /g s . However, the effects of such a superheavy 
DM will be eroded by inflation; moreover superheavy DM would not be produced significantly during a reheating 
phase of the Universe after its exit from the inflationary period. For intermediate string scales [27j , where the quantity 
M s /g s could be of order 10 11 GeV (which is the order of the GZK cutoff of ultra-high energy cosmic rays), there 
could be significant modifications in the relic abundances of superheavy DM particles with masses of this order. 
Such super-heavy DM particles can be produced during reheating [28^ . but in view of our scenario above, their relic 
abundance will be modified from the standard cosmology result. The presence of super-heavy DM, with increased relic 
abundances, might provide an explanation for the production of at least part of the spectrum of the ultra high energy 
cosmic rays, with energies of order 10 20 eV. Hence, the effects of D-foam on such scenarios are worthy of investigating 
further, especially in view of the fact that the density of D-particles might be significantly higher at earlier eras of 
the Universe, leading to stronger stochastic effects O(Af). 

For low string scales, of order a few TeV, the effects of the D-foam on thermal relic densities would be more 
significant. In fact, depending on the type of DM considered the effect could be constrained or falsified already by 
the WMAP five year data, since the induced increase of thermal relic abundance leaves less room for supersymmetry 
in the relevant parameter space. In certain cases, it may exceed the allowed region set by WMAP. The situation may 
thus lead to modifications of supersymmetry searches at colliders, especially in the context of neutralino DM models 
with the neutralino being Higgsino- or Wino-like, with masses up to TeV. In such cases there may be other reasons 
for an increased relic abundance, for instance slepton co-annihilation p9j . and in fact our effects of the foam are of 
comparable strength in some of these cases. 
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4 We remark at this stage that 11.261 has been derived by calculating the integral in 11.171 under the assumption that the scale factor 
a(t) evolves with temperature T as T ~ ^tj> wnere lb denotes today's temperature. However, this is only approximately true, since 
particles' annihilation at different stages of the evolution of the universe slowed down the cooling due to heat deposition. Actually, 
the corrected cooling law, obtained by assuming that the entropy (which is dominated by the relativistic species contribution) remains 

constant despite the presence of the source, reads: a(t) ~ ( s "//(*) ) a °T° ■ ^ ne f unc tion g e ff(t) is known [19] and is such that 
9e//,0 — 3-36 (due to the fact that today only photons and the light neutrinos contribute), while at freeze-out g e ff (tf) — 106.75 
within the framework of the Standard Model. In principle, one should calculate the integral in 11.17( numerically using this evolution 
of g B f f(t(T)). For our qualitative purposes in this work we shall ignore such corrections which, at any rate, do not affect significantly 
the order of magnitude of foam contributions to the relic abundances. To get an estimation of how big the difference is and assuming 
that the entropy will be conserved despite the presence of the source term (see also Il7ll ) implies that g e f fa 3 (t)T 3 will remain constant 
throughout the evolution of the universe (only relativistic species are connsidered to contribute to the entropy). Today only photons 
and neutrinos can contribute to the entropy yielding g e ff,o = 3.36. Therefore very roughly one would get a correction of the order of 

(a(t)T)| v ^ l Thig . g taken intQ accQunt ; n 
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